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Force correcting atom centered potentials for generalized gradient approximated 
density functional theory: Approaching hybrid functional accuracy for geometries and 
harmonic frequencies in small chlorofluorocarbons 
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^ Argonne Leadership Computing Faeility, Argonne National Laboratory, Argonne, Illinois 60439, USA 
^Department of Chemistry, University of Basel, Klingelbergstr. 80, 4056 Basel, Switzerland 

(Dated; March 22, 2013) 

Generalized gradient approximated (GGA) density functional theory (DFT) typically overesti- 
mates polarizability and bond-lengths, and underestimates force constants of covalent bonds. To 
overcome this problem we show that one can use empirical force correcting atom centered potentials 
(FCACPs), parameterized for every nuclear species. Parameters are obtained through minimiza- 
tion of a penalty functional that explicitly encodes hybrid DFT forces and static polarizabilities of 
reference molecules. For hydrogen, fluorine, chlorine, and carbon the respective reference molecules 
consist of H2, F2, CI2, and CH4. The transferability of this approach is assessed for harmonic 
frequencies in a small set of chlorofluorocarbon molecules. Numerical evidence, gathered for CF4, 
CCI4, CCI3F, CCI2F2, CCIF3, GIF, HF, HGl, CFH3, CF2H2, GF3H, CHGI3, GH2CI2, CH3CI indi- 
cates that the GGA-fFGACP level of theory yields harmonic frequencies that are significantly more 
consistent with hybrid DFT values, as well as slightly reduced molecular polarizability. 



I. INTRODUCTION 

The rigorous analysis of theoretical predictions and ex- 
perimental measurements for molecular properties is one 
of the noblest tasks in physical chemistry. The funda- 
mental importance of this line of research can hardly be 
overstated. To paraphrase M. Quack from one of his 
physical chemistry lectures at ETH Ziirich, we only un- 
derstand molecules when we are able to predict their prop- 
erties to a degree considered quantitative. Philosophically 
it is compelling to note that this definition distinguishes 
predictive power as a sufficient criterion for understand- 
ing. Conversely, qualitative predictions would indicate 
an only incomplete understanding. We also note that 
this definition qualifies "understanding" in terms of an 
arbitrarily chosen accuracy criterion for what, in some 
context of usefulness, is deemed sufficiently quantitative. 

Over the last decades, members of the laboratory of 
M. Quack have made major contributions to enable and 
carry out quantitative predictions for small and isolated 
molecules, using high-resolution molecular spectroscopy 
and highly accurate quantum chemical approaches to ad- 
dress fundamental questions such as parity violationii^. 
When it comes to condensed phase spectroscopy, how- 
ever, we ordinarily have to rely on less accurate gen- 
eralized gradient approximations (GGA) within density 
functional theory (DFT)^ii based molecular dynamics 
(AIMD)^. See for example the recent contributions by 
M. P. Gaigeot and others concerning infra red spectra 
of small peptides in liquid water |6|-|8|, as well as refer- 
ences therein. Vibrational spectra of surfaces can also be 
calculated using DFT—. And even for small molecules ad- 
sorbed on surfaces, the computational pediction of infra 
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red spectra has already been evinced, e.g. water adsorbed 
on Ni (111) and (211)iSiii, or organic molecules adsorbed 
on Si (111)1^. 

Unfortunately, GGAs are not always sufficiently accu- 
rate when it comes to vibrational properties, and admix- 
ture of Hartree-Fock exchange in the form of hybrid func- 
tionals might be necessary to yield significantly more ac- 
curate predictions^^, as also confirmed for the vibrational 
spectra of CCI3F— . For plane-wave basis calculations of 
condensed phase systems under periodic boundary con- 
ditions, however, the calculation of the exchange term is 
computationally dramatically more expensive than pure 
GGA, typically by an order of magnitude. One strategy 
to overcome this challenge consists of implementing so- 
phisticated and highly parallelized software that can effi- 
ciently exploit high-performance compute hardware such 
as IBM's BlucGenei^, with tens of thousands of compute 
nodes, e.g. CPMDi^iii or Qbojsi^iiS. While hybrid func- 
tionals are implemented and not impossible to use for 
condensed systems, as demonstrated in studies of liquid 
water—, and its IR spectrum^!, access to the substantial 
CPU resources required is as restricted as is the number 
of possible systems that can be tackled. 

At the hybrid DFT level of theory, comparative studies 
of a multitude of condensed systems are therefore pro- 
hibitive, let alone any exploration attempts in chemical 
compound space for bio or materials design from first 
principles'^. Alternatively, one can also attempt to im- 
prove the GGA's accuracy without increasing its compu- 
tational complexity. 

In this article, we will discuss and investigate a force 
correcting atom centered potentials (FCACPs) approach. 
The goal is to augment GGA calculations to approach 
hybrid DFT accuracy at negligible additional computa- 
tional cost. We will show that for all the molecules stud- 
ied, GGA-I-FCACP consistently yields geometries and 
harmonic frequencies with near hybrid DFT accuracy. 
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II. METHOD 

A. Background 

Within the Born-Oppenheimer approximation, 
HcUmann-Feynman's theorem clearly states that the 
forces on atoms are due to the electrostatic field exerted 
by the ground-state electron density^^. Since we ignore 
the explicit form of the exact exchange-correlation 
potential^ one can argue that it is reasonable to empiri- 
cally manipulate the density with the objective to yield 
forces that come as close as possible to known reference 
values. An appealing way to induce such changes into 
the electron density consists of adding atom centered 
non-local potentials that become negligibly small at the 
site of the nucleus, not to affect the core electrons, but 
are of sufficient magnitude in the interatomic region 
where the covalent bonding occurs. These correcting 
potentials are similar to pseudopotentials (PPs), and, 
while not necessarily so, can even assume their func- 
tional form. PPs, or effective core potentials, replace 
the explicit treatment of the core electrons^ir— , thereby 
reducing (i) the number of orbitals to be dealt with, and 
(ii) dramatically accelerating basis-set convergence in 
plane-wave basis sets. 

The idea to adapt PPs to account also for other prop- 
erties, i.e. going beyond the mere purpose of model- 
ing the core electrons' potential, is not new. It has 
successfully been deployed for relativistic effects^"—, 
self-interaction corrections modeling exact-exchange 
electron densities^, atomization energies and geome- 
tries of Al-clusters^, and minimizing quantum mechani- 
cal/molecular mechanical boundary error o'^^''^^" — , widen- 
ing the band gaf42r— , and introducing van-der-Waals 
interactions^^—. Various recent applications demon- 
strate the success of the latter, including binding of el- 
lipticinc to DNA and other biomolecular interactions^, 
and the accurate description of molecular crystals^i^ 
Furthermore, one can interpolate pseudopotentials, and 
perform self-consistent field calculations as a function 
of order parameter, < A < 1, effectively correspond- 
ing to fractional nuclear charges^"—. The effect of 
such "alchemical" variations on hydrogen-bonded dimcrs 
was investigated in combination with atom centered van 
der Waals correction^. Alchemical changes, and cor- 
responding Hellmann-Feynman derivatives^, are com- 
monly used for two, often related, purposes: Either for 
the evaluation of free energy differences between differ- 
ent compounds, e.g. using thermodynamic integration^, 
AF = J dX {dE/dX), e.g. sec Ref. |H; or for obtaining 
gradients that quantify a system's response to a variation 
in chemical compositio n^"^'^^!^"" — . 

B. Optimization 

The goal of this study is to explore if yet another 
property can be optimized through manipulation of 



PP parameter space, namely the force in the covalent 
bond. While at first GGAs can be considered suffi- 
ciently accurate, their interatomic distances, polarizabil- 
ities, and vibrational frequencies are typically consider- 
ably off when compared to more accurate methods, such 
as hybrid DFT results^^. While the choice of reference 
method and geometry is somewhat ad hoc, for this study 
the non-empirical hybrid functional PBEO has been se- 
lected^"— without any loss of fundamental generality: 
Post-Hartree-Fock methods such as Coupled Cluster, or 
quantum Monte Carlo, could have been chosen just as 
well. Homo-diatomics have been used as reference sys- 
tems, i.e. H2, F2, and CI2 for parameterizing the hy- 
drogen, fluorine, and chlorine atom. Since the smallest 
stable molecule consisting exclusively of carbon is the ex- 
ceedingly large C20 buckyball, methane has been chosen 
as a reference system instead. For the hydrogen atoms 
the previously optimized FCACPs are used without any 
further changes. 

For the optimization a unitless penalty functional, V , 
is defined that allows for multi-objective optimization of 
N normalized properties in PP parameters x, 

minT'fx) = min — V ' ' (1) 

where P^'^^ , P"^\ P^°^ are the instantaneous, initial, and 
reference values of property i, respectively. For this study 
Eq. ([1} has been chosen to minimize the deviation from 
two reference properties calculated with the hybrid func- 
tional PBEO^"— , (i) the 2-norm of all ionic forces of the 
system in the reference geometry, \F\, and (ii) the trace of 
the static polarizability tensor, Tr(a). The choice of the 
latter is based on the fact that the variational principle 
also holds in terms of Pearsson's maximum hardness prin- 
ciple^"— , an alternative to the potential energy whose 
origin is arbitrary in plane-wave PP based calculations. 

The minimization in PP parameter space could have 
been carried out with steepest descent or conjugate gra- 
dient algorithms using property derivatives calculated 
through linear response^S, as proposed in Ref. [s^. To 
facilitate the implementation, however, a gradient-free 
optimizer has been used for this study, namely Nelder- 
Mead's simplex optimization method^. The FCACP as- 
sumes the form of the highest empty angular momentum 
channel I for each atom type in the form of a Gocdecker- 
Hutter non-local PP— , in close analogy to the DCACP 
parameterization — , Specifically, I = s = 1 for H, and 
/ = c? = 3 for C, F, and CI. As such, this defines a 
2-dimensional parameter space of atom centered poten- 
tial parameters, x = (a;i,a;2), for each atom type. Here, 
xi — ri, the Gaussian width in the non-local projector 

plm 72 

{r\pi^) = NiYi^ifye'^, (2) 

where Ni and Yi,ri{f) correspond to the normalization 
constant and spherical harmonics, and r is the radial dis- 
tance from the atom. While this projector is centered on 
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TABLE I: Converged PBE+FCACP parameters for hy- 
drogen, fluorine, chlorine and carbon atoms, respectively 
optimized for reference systems H2, F2, CI2, and CH4 in 
PBEO geometry. |F|(P-|-F) denotes the residual Euclidean 
norm of PBE+FCACP forces in PBEO relaxed geometry in 
Hartree/A. A shows the deviation of trace of static polariz- 
ability from PBEO, for PBE (P) and PBE+FCACP (P+F), 
respectively. [Bohr"^] are shown as well. I = 1, 3, 3, 3, for H, 
F, CI, and C, respectively. 



Atom 


n [a.u.] 


hi [a.u.] 


|F|(P+F) A(P+F) A(P) 


H 


0.9871 


-0.004129 


3.8x10"' 


1.01 


1.15 


F 


1.3343 


-0.014713 


5.5x10"^ 


1.85 


2.02 


CI 


1.3199 


-0.006414 


3.9x10"^ 


5.94 


6.12 


c 


0.8081 


-0.035297 


3.8x10""* 


3.06 


3.43 



the atom, the multiplication of the polynomial term 
scales it down to zero at the position of the nucleus, xi 
thus tunes the location of the projector's maximum at a 
distance from the atom. The second parameter, X2 = hi, 
scales the magnitude of the entire nonlocal pseudopoten- 
tial contribution from the correcting channel. 



TABLE II: Results for four molecules used as reference sys- 
tems for optimization. Dipole moments, polarizabilities and 
harmonic frequencies with PBE, PBEO, and PBE+FCACP. 
PBE+FCACP results for PBEO geometries. Traces of polariz- 
abilities |a| in atomic units. Harmonic frequencies uj [cm~*] in 
descending order of energy. If degenerate the averaged value 
is reported, d and t denoting doublet or triplet. Experimental 
wavenumbers u, all in [cm~*] and from Ref.— , are given as a 
footnote for orientation. 



Molecule 


method 


|a| 








CJ4 


Ha" 


PBE 
PBE+FCACP 
PBEO 


6.84 
6.66 
4.40 


4306.6 
4386.2 
4402.7 










PBE 
PBE+FCACP 
PBEO 


12.99 
12.49 
6.80 


1087.2 
1155.2 
1139.2 










PBE 
PBE+FCACP 
PBEO 


41.92 
41.47 
24.56 


538.6 
561.2 
578.1 








CH4'' 


PBE 


16.87 


3075.56* 


2951.3 


1494.3" 


1277.5* 




PBE+FCACP 


16.50 


3142.0* 


3017.9 


1514. 7^* 


1287.8* 




PBEO 


13.44 


3157.2* 


3029.8 


1539.7^* 


1316.9* 



"4401.2; 916.9; 559.7; ''SOig.O, 2917, 1534.0, 1306.0 



III. RESULTS AND DISCUSSION 

VFCACp(i",r) = 2^{r\pim)hi{pi,n\r ). (3) 

Im 



C. Computational details 

All the DFT reference and optimization calculations 
have been carried out with the CPMD PP plane-wave ba- 
sis set prograro^^. Polarizabilities are obtained through 
the linear response tools implemented in CPMoSi, and har- 
monic frequencies in CPMD are calculated from Hessians 
obtained via finite differences. PBE Goedecker-Hutter 
PPs have been used for all the calculation o^^i^'^ . While 
the use of GGA PPs within hybrid DFT, or even within 
GGA + atom centered corrections, is not unconventional, 
we note that eventually the PPs should be reparame- 
terized to be entirely consistent with their density func- 
tional2^. All calculations involved isolated boundary con- 
ditions with the Poisson solver by Martyna and Tuck- 
erman2^, a plane-wave cutoff of 200 Ry, a unit-cell of 
15x15x15 A"^, and were carried out on Argonne Leader- 
ship Computing Facility's IBM BlueGene/P machine. At 
this point the reader is cautioned that severe finite size 
effects in the plane-wave calculations can lead to signifi- 
cant distortions of polarizabilities and frequencies. This 
is of little concern for this study, however, since we only 
deal with relative changes in these properties, and the 
finite size effects can be assumed to cancel when compar- 
ing results from different functionals. 



A. Optimization for reference molecules 

Optimized FCACP parameters are reported in Table U 
for hydrogen, fluorine, chlorine, and carbon. The po- 
sitioning parameter, rj, converges to a length scale of 
0.8 to 1.4 Bohr for all atoms. For comparison, the ac- 
tual PP parameters in the non-local PP s-channel (not 
existent for the hydrogen's PP) range from 0.2 to 0.3 
Bohr for C, F, and CI. By contrast, the DCACP values 
for the /-channel range from 1.8 to 3.6 Bohr— i^^. It is 
thus clear that the FCACP corrections acts in a more 
intermediate mid-range around the atoms. This range of 
action is not surprising, it is roughly the distance from 
atom to covalent bonding. Similar comparisons can be 
made for the magnitude of the correction, hi: The s- 
channel of the C, F, and CI PPs ranges from 9.6 to 23.7 
a.u. in magnitude, the DCACP's /-channel is typically 
of only - -10-4 a.uM^. The FCACP's hi is in between 
at ^ -10"^ to -lO"'^ a.u. It is not surprising that the 
FCACPs are one to two orders of magnitude larger than 
the DCACPs since the absolute errors in covalent forces 
are significantly larger than the GGA's error to account 
for van der Waals forces. As such, FCACPs correspond 
to an attractive potential that counteracts the effect of an 
overly delocalized GGA electron density on interatomic 
covalent forces. 

Table |T] also enlists the residual forces and polariz- 
abilities in the corresponding reference molecules after 
convergence of the FCACPs. We note that on the one 
hand the norm of the forces can be quenched to virtu- 
ally correspond to the reference method's force (zero in 
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this study) within numerical precision, with the shght 
exception of carbon where convergence sets in aheady at 
-10"^ Hartree/A. The uncorrected PBE forces in PBEO 
geometry typically amount to 10~^ to 10"'^ Hartree/A. 
On the other hand, however, the static polarizability's 
trace does not improve as much, the PBE value is usu- 
ally reduce by no more than 10%. Additional optimiza- 
tion test runs suggest that removal of the force from the 
objective penalty functional in Eq. ([T]) would yield per- 
fect reproduction of the reference polarizability, at the 
expense, however, that the resulting forces worsen con- 
siderably. As such, polarizability appears to be crucial 
as a constraint for the main objective of having correct 
forces in the covalent bond. This finding underscores the 
importance of polarizability, as already widely researched 
and discussed in terms of the maximal hardness princi- 
ple and the hard-soft acid-base principles by Pearson, 
and Parr and Chattaraj^"— . 



B. Frequencies for reference molecules 

Tabic |TT] reports final polarizability traces, as well as 
harmonic frequencies for the four reference molecules 
used for optimization, i.e. H2, F2, CI2, and CH4. We 
note that even though frequencies were not explic- 
itly encoded in the penalty minimization, for all cases 
the PBE-I-FCACP frequencies are significantly closer 
to PBEO. For all frequencies except the one for F2, 
PBE-fFCACP approaches the PBEO value from below. 
Already this suggests that F containing molecules may 
be harder to deal with, vide infra. 



C. Transferability to other molecules 

The transferability of the PBE-fFC AGP functional has 
been tested for 14 other molecules and properties. Specif- 
ically, Table IIIII gives an overview for all the calculated 
dipole moments, polarizabilities, and harmonic frequen- 
cies for all test molecules. As one would expect, the hy- 
brid functional PBEO yields consistently higher frequen- 
cies, lower polarizabilities, and lower dipole moments 
than the GGA functional PBEii^. For most molecules, 
the FGAGP corrected GGA results indicate a clear and 
systematic shift towards the hybrid functional harmonic 
frequency numbers; not only for the molecules used for 
training (Table HI]) but also for all the "unseen" test 
molecules: The mean absolute deviation from the PBEO 
results for the test molecules reduces dramatically, in 
case of the highest lying mode even from 100 to 33 
cm~^. As the frequency of the mode decreases, the de- 
viation from PBEO becomes smaller for both, PBE and 
PBE-I-FCACP. Eventually, for the lowest lying modes, 
the correction performs slightly worse than PBE, but 
is still within reason. Encouragingly, we also note that 
PBE+FCACP rarely overestimates the PBEO frequency, 
the largest overcstimation being ^28 cm~^ in the case of 



in CCI2F2. 

As already alluded to above, the correction performs 
worst when Fluorine is present, albeit not always. For ex- 
ample, LOi of CH2F2 is modeled worse by PBE+FCACP 
than by PBE alone, this is possibly also due to the fact 
that the PBEO frequency is already smaller than the PBE 
frequency in this case. But also for CF4 results for modes 
L02, ws, and 0-14 suggest that there is still room for im- 
provement for the F correction. 

For polarizabilities, the improvement is significantly 
less dramatic. As one would already expect from the op- 
timization penalty results discussed in the previous sec- 
tion (Table m, the diagonal polarizability tensor elements 
change only in the single digit percentages. The change, 
however, is consistently towards a reduction in polariz- 
ability, as desired. The PBE+FCACP dipole moments 
worsen when compared to PBE, albeit ever so slightly. 

It remains to be seen if PBE+FCACP can also ac- 
count for anharmonicities. In analogy to using multiple 
DCACP channels to generate the correct dissocia- 
tive behavior in the case of rare gas atoms^, additional 
FGAGP channels might be necessary. The PBE+FCACP 
forces on all atoms placed in the PBEO geometry are typ- 
ically rather small, but not negligible. No imaginary fre- 
quencies have been found for any of the molecules investi- 
gated. A comprehensive assessment of these and further 
properties, however, is beyond the scope of this study, 
and will be the topic of future work. 



IV. CONCLUSION 

A force correcting atom centered potential (FGAGP) 
has been introduced for augmenting generalized gradient 
approximated (GGA) DFT calculations. FGAGPs have 
been trained and tested using hybrid DFT reference cal- 
culations of various small molecules involving hydrogen, 
chlorine, fluorine and carbon atoms. The optimization of 
parameters has been accomplished by minimizing forces 
on atoms in hybrid DFT relaxed geometries, as well as de- 
viation of static polarizabilities from hybrid DFT values. 
Numerical evidence suggests that not only hybrid DFT 
geometries can be achieved, but that also the correspond- 
ing harmonic frequencies improve dramatically when us- 
ing GGA+FGAGP. The promising numerical results and 
the versatility of approach hold great promise that vi- 
brational studies of condensed phase systems are possible 
with hybrid DFT accuracy — at the computational cost 
of GGA calculations. 

In the future, it might be worthwhile to more system- 
atically investigate the construction of better functionals 
using series expansions of atom centered potentials on 
top of typical GGAs^^, and to assess the effect on multi- 
ple properties at once. For example, one could combine 
FGAGPs with DCACPs to calculate vibrational prop- 
erties in molecular liquids or crystals, as well as phase 
diagrams. It remains to be seen which of the many ap- 
proximations to the exchange-correlation potentials are 
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TABLE III: Results for test molecules for which the FCACPs have been used without any further changes. Dipole moments, 
polarizabilities and harmonic frequencies with PBE, PBEO, and PBE+FCACP. PBE+FCACP results for PBEO geometries. 
Total dipole moments /i and averaged traces of polarizability |a| in atomic units. Harmonic frequencies ui [cm~^] in descending 
order of energy. If degenerate the averaged value is reported, d and t denoting doublet or triplet. Experimental wavenumbers 
u, all in [cm~^] and from Ref. [ttI . are given as a footnote for orientation. The mean absolute error from PBEO using either 
PBE or PBE+FCACP is shown in the last two lines, respectively. 



Molecule 


method 




\a\ 




UJ2 


W3 


LtJ4 




We 








HF " 


PBE 

PBE+FCACP 
PBEO 


0.71 
0.68 
0.72 


6.02 
5.99 
4.69 


3932.2 
4233.2 
4231.9 


HCl' 


PBE 
PBE+FCACP 
PBEO 


0.43 
0.41 
0.44 


17.58 
17.52 
14.17 


2900.4 
2989.0 
2999.5 


CIF'^ 


PBE 

PBE+FCACP 
PBEO 


0.33 
0.31 
0.33 


18.64 
18.47 
14.75 


791.0 
840.6 
852.7 


ecu'' 


PBE 
PBE+FCACP 
PBEO 


0.02 
0.02 
0.02 


70.40 
69.85 
55.46 


729.1* 
745.6* 
799.4* 


440.1 
445.4 
467.6 


300.4 
306.8* 
315.6' 


199.6" 
204.5" 
212.1" 












CF4" 


PBE 
PBE+FCACP 
PBEO 


0.01 
0.01 
0.01 


20.38 
20.04 
15.91 


1230.1* 
1291.6* 


9UD. I 

862.7 
930.9 


DO/ .U 

639.5* 
677.9* 


Odd. 9 

496.4" 

539.6" 












CClsF^ 


PBE 
PBE+FCACP 
PBEO 


0.19 
0.19 
0.19 


57.83 
57.34 
45.47 


1046.2 
1053.4 
1124.5 


797.5" 
826.8" 
862.1" 


537.1 
541.4 
563.6 


458.1" 
464.6" 
463.4" 


354.8 
359.8 
369.5 


250.0" 
257.3" 
262.2" 








CHCl3« 


PBE 
PBE+FCACP 
PBEO 


0.40 
0.37 
0.42 


57.11 
56.66 
45.02 


3066.3 
3138.7 
3159.9 


1187.4" 
1200.7" 
1252.5" 


722.3" 
734.9" 
788.9" 


657.2 
670.3 
689.7 


352.7 
360.5 
371.8 


245.4" 
250.4" 
257.9" 








CFH3" 


PBE 
PBE+FCACP 
PBEO 


0.73 
0.70 
0.73 


17.48 
17.14 
13.74 


3055.5" 
3133.8" 
3143.8" 


2946.1 
3023.9 
3034.8 


1531.8" 
1556.8" 
1553.4" 


1481.9 
1494.7 
1522.6 


1291.9" 
1331.5" 
1318.1" 


1065.9 
1062.8 
1126.8 








CCIF3' 


PBE 
PBE+FCACP 
PBEO 


0.19 
0.18 
0.19 


32.60 
32.17 
25.54 


1147.7" 
1169.0" 
1231.3" 


1051.5 
1070.8 
1114.7 


799.0 
770.3 
800.3 


639.3" 
568.0" 
584.3" 


483.3 
455.5 
480.5 


446.8" 
355.1" 
365.0" 








CHsCl-' 


PBE 
PBE+FCACP 
PBEO 


0.74 
0.71 
0.75 


29.67 
29.33 
23.65 


3084.8" 
3154.2" 
3174.2" 


2983.2 
3053.0 
3069.0 


1419.2" 
1434.2" 
1468.3" 


1326.3 
1336.7 
1376.7 


999.1" 
1014.4" 
1038.7" 


716.7 
720.7 
757.3 








CHFa" 


PBE 
PBE+FCACP 
PBEO 


0.62 
0.58 
0.64 


19.63 
19.29 
15.25 


3074.5 
3153.9 
3158.3 


1527.7" 
1524.3" 
1522.3" 


1139.2" 
1153.5" 
1212.4" 


1122.9 
1133.6 
1182.5 


755.5 
687.4 
762.0 


628.7" 
501.7" 
613.9" 








CCI2F2' 


PBE 
PBE+FCACP 
PBEO 


0.21 
0.21 
0.22 


45.13 
44.69 
35.44 


1102.1 
1119.3 
1189.4 


1058.0 
1076.6 
1133.5 


867.0 
906.6 
931.3 


723.7 
703.7 
714.7 


499.6 
500.3 
500.5 


492.7 
484.1 
492.2 


460.8 
469.3 
479.3 


440.0 275.9 
459.7 278.5 

432.1 280.6 


CHaCb™ 


PBE 
PBE+FCACP 
PBEO 


0.62 
0.59 
0.64 


43.37 
42.99 
34.33 


3099.6 
3169.5 
3190.3 


3021.7 
3093.0 
3111.5 


1401.4 
1416.9 
1454.2 


1226.7 
1239.9 
1283.0 


1139.0 
1154.7 
1186.1 


876.3 
897.1 
912.2 


720.9 
728.6 
783.4 


693.6 
701.6 
732.3 


263.5 
274.3 
280.9 


CH2F2" 


PBE 
PBE+FCACP 
PBEO 


0.76 
0.72 
0.77 


18.50 
18.19 
14.42 


3073.0 
3113.4 
3159.8 


2975.8 
3043.2 
3066.4 


1573.7 
1495.7 
1598.0 


1554.5 
1419.7 
1547.4 


1524.5 
1240.1 
1522.3 


1271.7 
1192.3 
1288.3 


1133.0 
1080.7 
1186.5 


1078.5 
1043.0 
1155.5 


678.8 
567.1 
661.9 


MAE 
MAE 


PBE 
PBE+FCACP 


0.01 
0.03 


6.89 
6.56 


100.1 
33.3 


61.9 
31.7 


36.7 
37.6 


32.3 
35.1 


17.7 
52.5 


30.6 
39.4 


44.8 
56.9 


41.2 
56.9 


13.0 
34.5 



"4138.4; ''2990.9;<=738.5; '^776.0, 459.0, 314.0, 217.0; =1280.0, 
909.0, 631.0, 453.0; ■''1085.0, 847.0, 535.0, 394.0, 350.0, 241.0; 
93034.1, 1219.7, 774.0, 680.0, 366.0, 260.0; ''3006.0, 2930.0, 1467.0, 
1464.0, 1182.0, 1049.0; ^212.0, 1105.0, 781.0, 563.0, 476.0, 350.0; 
^3041. 8, 2966.2, 1454.6, 1354.9, 1015.0, 732.1; '=3036.0, 1372.0, 
1152.0, 1117.0, 700.0, 507.0; '1159.0, 1101.0, 902.0, 667.0, 458.0, 
446.0, 437.0, 322.0, 262.0; ™3040.0, 2999.0, 1467.0, 1268.0, 1153.0, 
898.0, 758.0, 717.0, 282.0, "3014.0, 2948.0, 1508.0, 1435.0, 1262.0, 
1178.3, 1111.2, 1090.1, 528.5 
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the most suited for being combined with atom centered 
potentials. Apart from PBE, the authors also consid- 
ers the AMOS functional^ to be a potential candidate 
for such an effort. We do not know yet if such an 
GGA+FCACP+DCACP approach can account for all 
the many-body interactions present, recently found to 
be relevant even for van der Waals interactions^^Ti^. An- 
other critical issue is to more rigorously address the some- 
what arbitrary choice of reference systems (molecules and 
geometries) and methods (other than hybrid DFT). Re- 
cent efforts using machine learning in chemical compound 
space might offer new strategies to remove the resulting 
selection bias^ 
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